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1.  Introduction 


As  the  focus  of  technology  continues  to  become  smaller,  the  need  for  tools  to  accurately  image 
nanoscale  samples  continues  to  grow.  A  large  step  toward  meeting  this  need  occurred  in  the 
1980s  with  the  development  of  atomic  force  microscopy  (AFM).  With  the  invention  of  AFM,  the 
ability  to  image  samples  on  the  atomic  scale  became  more  accurate  and  reliable.  Over  the  past  20 
years,  numerous  changes  have  been  made  in  the  different  modes  of  operation  and  materials  used 
during  scanning.  These  changes  have  led  to  increased  performance  in  atomic  scale  imaging  and 
have  consistently  increased  the  potential  applications  of  AFM.  Along  with  determining  surface 
topography  of  samples,  AFM  has  been  used  to  obtain  information  about  the  material  and  chemical 
properties  of  samples.  The  versatility  in  type  of  information  gathered  through  AFM  has  made  it 
one  of  the  most  effective  tools  in  imaging  on  the  atomic  and  nanoscale  levels. 

The  atomic  force  microscope  (see  figure  1)  is  centered  around  a  microscale  cantilever  beam  that 
has  a  shape  tip  at  the  end.  These  tips  are  usually  made  of  silicon  or  carbon  nanotubes,  depending 
on  the  desired  application.  Although  each  mode  of  operation  varies,  the  basic  principles  are  the 
same.  The  cantilever  beam  is  excited  at  a  certain  amplitude  and/or  frequency.  The  interaction 
forces  between  the  tip  and  the  sample  lead  to  changes  in  the  amplitude,  phase,  or  frequency  of 
oscillation.  Information  can  be  obtained  about  the  sample  because  of  these  changes.  The  motion 
of  the  cantilever  beam,  along  with  the  changes  in  motion,  is  often  determined  by  a  diode  laser 
reflected  onto  an  array  of  photodiodes.  The  motion  of  the  tip  can  be  fit  to  a  sine  curve  that  informs 
the  control  loop  of  the  amplitude,  frequency,  and  phase  of  the  cantilever  at  all  times. 


Figure  1.  Basic  operation  of  an  atomic  force 

microscope  (http://www.physics.ucsb. 
edu/~hhansma/afm-acs_news.htm). 
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The  governing  equation  for  the  vibration  of  the  beam  is  given  by 


m  Z+  kZ  + 


mco0 

~Q 


z  = 


F,  +  Fo  cos  (cat) . 


(1) 


This  equation  considers  the  cantilever-probe  system  as  a  point-mass  spring  system  in  which  m  is 
the  mass  of  the  cantilever,  k  is  the  cantilever  force  constant,  co0  is  the  free  resonance  frequency,  Q 
is  the  quality  factor,  Fts  is  the  tip-sample  interaction  forces,  F0  is  the  amplitude  of  the  excitation 
force,  co  is  the  frequency  of  the  excitation  force,  Z  is  the  height  of  the  probe,  and  t  is  time. 


2.  Tip  Sample  Interaction  Forces 


As  the  tip  or  probe  moves  relative  to  the  sample,  it  experiences  varying  forces.  These  forces 
include  but  are  not  limited  to  Van  der  Waals  forces,  mechanical  contact  forces,  chemical  bonding, 
electrostatic  forces,  and  capillary  forces.  The  tip-sample  interaction  forces  depend  on  the  type  of 
tip  used,  the  material  of  the  type,  and  the  material  of  the  sample.  These  forces  can  be  modeled  in 
a  number  of  ways  such  as  with  the  Lennard- Jones  potential,  which  is  a  relatively  simple  model 
representing  the  interaction  between  molecules.  The  force  attributable  to  the  Lennard- Jones 
potential,  which  is  the  negative  gradient  of  the  potential  energy,  is  represented  by 


F(r)  =  4  s 


f  CT12  CJ6A 
12  — -6  — 

V  r  r  j 


(2) 


where  r  is  the  distance  between  molecules,  s  is  the  depth  of  the  potential  well,  and  a  is  the 
location  of  zero  potential. 

Figure  2  shows  the  shape  of  a  Lennard-Jones  force  curve.  When  the  tip  and  sample  are  far  away 
from  each  other,  there  is  a  small  attractive  force.  As  the  tip  moves  downward  toward  the  sample, 
the  force  becomes  more  attractive  until  it  reaches  the  maximum.  As  the  tip  continues  to  move 
down,  the  force  becomes  less  attractive  until  it  reaches  the  repulsive  regime. 


For  the  simulations  in  this  work,  the  tip-sample  interaction  curves  are  detennined  through  the 
calculation  of  the  discrete  atomic  interactions.  The  energy  of  the  system  is  calculated  as  a  probe  is 
lowered  toward  the  sample  at  a  fixed  horizontal  position.  We  find  the  force  by  taking  the  negative 
gradient  of  the  energy,  which  is  a  function  of  the  probe’s  distance  from  the  sample.  This  force  is 
then  fit  to  equation  3  for  use  in  the  frequency  and  force  modulation- ATM  (FFM-AFM)  simulation. 


F 


-W 


1  +  30  *(Z„-a) 

F  =  -W  +  S*(Z,,-<t)2 


for  Zts  >  a 
for  Zts  <  cr 


(3) 
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(a)  15-nm-radius  silicon  tip  approaching  a  Si(100)-OH  surface,  (b)  a  5.4-nm-radius  single-walled  carbon 
nanotube  tip  approaching  a  Si(100)-OH  surface,  and  (c)  a  2.4-nm-radius  triple-walled  carbon  nanotube 
approaching  a  bacteriorhodopsin  molecule 


Figure  2.  Illustrating  the  general  shape  of  a  force  versus  distance  graph  with  the  Lennard-Jones  potential  model  (7). 

In  equation  3,  Zts  represents  the  distance  between  the  tip  and  sample,  W  is  the  maximum  attractive 
force,  a  is  the  location  at  which  the  maximum  attractive  force  occurs,  and  S  is  the  slope  of  the 
curve  in  the  repulsive  regime. 

There  are  types  of  cantilever  beam  modulation  approaches  in  AFM.  One  mode  is  amplitude 
modulation  (AM-AFM)  which  is  also  referred  to  as  tapping  mode.  In  this  mode,  the  cantilever  is 
excited  at  or  near  resonance  frequency  at  a  height  well  above  the  sample.  The  cantilever  is  then 
lowered  toward  the  sample  while  the  excitation  frequency  and  amplitude  are  kept  constant.  As 
the  tip  enters  the  repulsive  regime,  the  effective  amplitude  of  the  tip  decreases.  The  height  of  the 
sample  is  determined  when  the  amplitude  reaches  a  set  point,  which  is  defined  before  the  experi¬ 
ment  and  is  below  the  freely  oscillating  amplitude.  There  are  two  main  weaknesses  with  this 
mode.  The  first  is  that  it  usually  results  in  relatively  high  tip-sample  interaction  forces.  This  can 
lead  to  damage  of  the  sample  if  the  sample  is  too  soft.  The  second  is  that  it  can  also  be  subject  to 
bistability  which  leads  to  inaccurate  height  measurements.  Bistability  occurs  when  the  tip  jumps 
from  the  attractive  regime  to  the  repulsive  regime,  and  vice  versa,  as  illustrated  in  figure  3. 


Figure  3.  Illustration  of  bistability  effect  in  AM-AFM. 
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The  second  common  mode  of  imaging  in  AFM  is  frequency  modulation- AFM  (FM-AFM).  In  this 
mode,  the  excitation  amplitude  is  kept  constant  but  the  excitation  frequency  is  varied  to  match  the 
current  resonance  frequency.  This  frequency  changes  because  of  the  tip-sample  interaction  by 


/  = 


k-VF„ 


m 


(4) 


where/is  the  current  resonance  frequency,  k  is  the  force  constant,  m  is  the  mass  of  the  cantilever, 
and  VFts  is  the  gradient  of  the  tip-sample  interaction  forces  over  the  range  in  which  the  tip 

is  oscillating. 

In  FM-AFM,  the  tip  height  is  decreased  until  the  desired  frequency  shift  is  attained.  The  tip  is  then 
moved  along  the  sample  and  continues  to  find  the  height  that  achieves  the  desired  frequency  shift. 
The  amplitude  of  oscillation  is  much  smaller  than  in  AM- AFM.  The  weakness  in  this  mode  of 
AFM  is  that  the  small  oscillations  can  have  difficulties  in  imaging  sharp  topographical  changes 
and  the  tip  does  not  get  as  close  to  the  sample,  which  may  mean  that  a  true  scan  is  not  achieved. 
This  effect  is  often  referred  to  as  tip  broadening  and  is  shown  in  figure  4. 


Figure  4.  Tip  broadening  because  the  tip  never  reaches  the  true  sample  skin. 


3.  Frequency  and  Force  Modulation-Atomic  Force  Microscopy  (FFM-AFM) 


Although  multiple  modes  of  AFM  have  been  developed  and  numerous  adaptations  made,  AFM 
still  struggles  to  accurately  image  soft  samples  without  causing  permanent  damage  to  the  sample. 
In  this  proposed  novel  mode  of  AFM,  the  cantilever  equilibrium  position  is  kept  constant.  The 
excitation  frequency  is  that  of  the  current  resonance  frequency,  as  in  FM-AFM.  The  excitation 
amplitude  is  varied  until  the  frequency  shift  from  the  free  resonance  frequency  is  zero.  At  this 
point,  the  attractive  and  repulsive  forces  are  equal.  The  sample  height  is  then  detennined  as  the 
cantilever  equilibrium  minus  the  amplitude  of  oscillation.  The  cantilever  sweeps  horizontally 
across  the  sample,  providing  a  real-time  scan  of  the  sample  surface.  The  control  scheme  for  FFM- 
AFM  is  depicted  in  figure  5. 
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Figure  5.  Control  scheme  for  FFM-AFM. 

FFM-AFM  has  many  benefits  over  other  AFM  modes,  especially  in  the  scanning  of  soft  samples 
such  as  biological  materials.  One  advantage  is  that  FFM-AFM  limits  the  maximum  repulsive 
force  because  of  tip-sample  interaction  forces  by  limiting  the  frequency  to  no  higher  than  the  free 
resonance  frequency.  Figure  6  illustrates  how  the  maximum  repulsive  force  using  FFM-AFM  is 
less  than  that  of  AM-AFM.  Another  benefit  is  that  bistability  is  eliminated  because  the  excitation 
frequency  is  always  the  same  as  the  current  resonance  frequency.  A  third  advantage  is  that  the 
control  scheme  includes  only  one  control  loop,  meaning  it  should  be  a  relatively  simple  method  to 
implement. 


Figure  6.  Simulation  of  a  triple-walled  carbon  nanotube  imaging 
a  bacteriorhodopsin  molecule  on  a  Si(100)-OH  surface. 
(The  square  indicates  the  maximum  repulsive  force 
experienced  during  AM-AFM  imaging.  The  circle  is 
the  force  experienced  during  FFM-AFM.) 
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4.  Quantitative  Prediction  Tip  Penetration  and  Repulsive  Forces 


The  goal  of  this  project  is  to  quantitatively  predict  the  maximum  tip-sample  repulsive  forces  and  the 
maximum  tip  penetration,  based  on  a  proposed  set  of  AFM  scanning  parameters.  These  parameters 
are  the  cantilever  force  constant  k,  the  quality  factor  Q,  the  fixed  equilibrium  position  of  the  canti¬ 
lever  Zc,  and  the  free  resonance  frequency  w0.  The  free  resonance  frequency  and  force  constant 
have  an  approximately  linear  relationship  and  therefore  were  varied  together.  We  performed 
numerical  simulations  to  model  the  imaging  of  a  single  scan  point  by  performing  numerical  inte¬ 
gration  of  the  equation  of  motion  (equation  1)  using  the  Verlet  algorithm  to  detennine  the  velocity 
and  position  at  each  time  step.  After  the  beam  achieves  a  zero  frequency  shift  from  the  free  reso¬ 
nance  frequency,  the  maximum  tip  penetration  is  determined  by  the  difference  between  the  lowest 
point  reached  by  the  probe  and  the  true  height  of  the  sample.  The  maximum  interaction  force  can  be 
calculated  when  this  value  is  entered  into  equation  3. 

The  base  parameters  used  for  each  run  are  listed  in  table  1 .  Table  2  shows  the  range  used  for  each 
parameter  and  the  resulting  variation  in  tip-sample  repulsive  force  and  tip  penetration  as  each  AFM 
scanning  parameter  is  varied. 


Table  1.  Base  parameters  used  to  investigate  the  correlation  between  AFM  scanning 
parameters  and  tip  penetration  and  repulsive  forces. 


Parameter 

Imaging  in  air 

Imaging  in  water 

Cantilever  force  constant,  N/m 

2.5 

2.5 

Cantilever  free  resonant  frequency,  kHz 

100 

33 

Cantilever  quality  factor 

150 

3 

Tip  rest  position  above  the  substrate,  nm 

10 

10 

Well  depth  (for  tip-sample  interaction  forces) 

0.5 

0.5 

Steepness  (for  tip-sample  interaction  forces) 

5 

5 

Table  2.  Results  from  numerical  simulations. 


Parameter 

Range 

Maximum  force  over  parameter 
range  (nano-newton  [nN]) 

Minimum  force  over 
parameter  range  (nN) 

Variation 

(percent) 

Air 

Liquid 

Air 

Liquid 

Air 

Liquid 

Zc,  nano¬ 
meters  (nm) 

5  to  50 

0.658 

0.652 

0.637 

0.638 

3.3 

2.2 

k,  newtons  per 
meter  (N/m) 

0.5  to  50 

0.669 

0.670 

0.639 

0.598 

4.8 

11.9 

Q 

5  to  200 

0.648 

0.648 

0.646 

0.646 

0.3 

0.4 

As  seen  in  table  2,  the  AFM  scanning  parameters  only  have  a  small  effect,  if  any,  on  the  maximum 
tip-sample  repulsive  forces.  The  variation  in  the  force  constant  showed  the  only  real  change  and 
this  resulted  in  a  0.07  nN  difference  in  force.  All  forces  are  well  below  the  critical  level  of  2.4  nN 
which  would  result  in  pennanent  damage  to  the  sample  (7).  The  variation  on  the  tip  penetration 
was  negligible  and  remained  between  0.47  and  0.48  ntn. 
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The  results  of  the  numerical  simulations  suggest  that  these  scanning  parameters  do  not  greatly 
affect  the  repulsive  forces  and  tip  penetration.  This  claim  can  be  justified  analytically  if  we  exam¬ 
ine  the  frequency  shift  equation  from  perturbation  theory  (2).  This  equation  assumes  the  cantilever 
acts  as  a  weakly  perturbed  hannonic  oscillator.  Although  this  assumption  is  only  a  true  representa¬ 
tion  of  the  oscillating  cantilever  beam  when  the  restoring  force  of  the  beam  is  significantly  greater 
than  the  interaction  forces,  it  is  a  good  basis  for  analytical  study.  When  we  define  d+A0  as  Zc  and 
set  the  frequency  shift  to  zero  in  equation  5,  the  equation  simplifies  to  equation  6. 


\  co 

A/  =  jf-  J Fts[d  +  A0+A0  cos((p)]cos((p)d(p 
Ln  kA0  0 

2k 

0  =  J Fts[Zc  +  A0  cos(»]cos {(p)d(p 

o 


(5) 

(6) 


As  illustrated  in  equation  6,  the  frequency  shift  is  not  a  function  of  the  quality  factor  (Q).  Also,  the 
force  constant  (k)  and  the  free  resonance  frequency  (co0)  are  removed  because  the  frequency  shift  is 
zero.  Therefore,  it  can  be  concluded  that  none  of  those  three  factors  have  an  impact  on  the  frequen¬ 
cy  shift  and  consequently,  the  maximum  interaction  forces.  By  performing  numerical  integration 
using  the  mid-point  method  with  10,000  intervals,  while  varying  Zc,  we  can  find  the  repulsive 
force.  With  this  infonnation,  the  tip  penetration  can  also  be  detennined.  Table  3  contains  the 
results  of  this  integration.  Both  the  repulsive  force  and  tip  penetration  remain  relatively  small  and 
unchanged  except  at  an  equilibrium  height  of  4  nm.  The  differences  between  these  data  and  the 
results  of  the  numerical  simulations  can  be  attributed  to  the  fact  that  equation  5  assumes  a  weakly 
perturbed  harmonic  oscillator. 


Table  3.  Results  from  numerical  integration  of  equation  6. 


W  depth 

0.5 

W  location 

3.4 

S 

5 

Zc 

Ao 

Zmin 

Fts(Zmin) 

Penetration 

4 

1.0240 

2.9760 

0.3991 

0.4240 

5 

2.0403 

2.9597 

0.4695 

0.4403 

6 

3.0462 

2.9538 

0.4953 

0.4462 

7 

4.0493 

2.9507 

0.5092 

0.4493 

8 

5.0513 

2.9487 

0.5182 

0.4513 

9 

6.0527 

2.9473 

0.5245 

0.4527 

10 

7.0537 

2.9463 

0.5293 

0.4537 

11 

8.0545 

2.9455 

0.5330 

0.4545 

12 

9.0552 

2.9448 

0.5360 

0.4552 

13 

10.0557 

2.9443 

0.5385 

0.4557 

14 

11.0562 

2.9438 

0.5406 

0.4562 

15 

12.0566 

2.9434 

0.5424 

0.4566 

16 

13.0569 

2.9431 

0.5440 

0.4569 

17 

14.0573 

2.9427 

0.5454 

0.4573 

18 

15.0575 

2.9425 

0.5466 

0.4575 

19 

16.0578 

2.9422 

0.5478 

0.4578 

20 

17.0580 

2.9420 

0.5487 

0.4580 

1.15356 

37.50489 

8.0034 

Percent  Change 

0.59948 

16.87206 

4.0052 

Percent  Change 
for  Zc  =  5-20 
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5.  Conclusions 


AFM  is  a  powerful  tool  with  numerous  applications  for  the  imaging  and  analysis  of  nanoscale 
samples.  Despite  the  wide  range  of  AFM  methods,  the  imaging  of  soft  samples  is  still  limited. 
FFM-AFM  is  a  mode  of  AFM  that  allows  for  the  accurate  scanning  of  soft  samples  without  the 
risk  of  sample  damage.  FFM-AFM  appears  to  be  nearly  independent  of  the  AFM  scanning 
parameters.  More  work  is  necessary  to  examine  more  parameter  combinations  and  to  begin 
determining  more  quantitative  ways  of  comparing  the  AFM  parameters  to  the  tip  penetration  and 
tip-sample  repulsive  forces.  A  further  investigation  into  other  analytical  methods  to  detennine  the 
correlation  is  also  required. 
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